Atomic scale lattice distortions and domain wall profiles 
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We present an atomic scale theory of lattice distortions using strain related variables and their 
constraint equations. Our approach connects constrained atomic length scale variations to contin- 
uum elasticity and describes elasticity at several length scales. We apply the approach to a two- 
dimensional square lattice with a monatomic basis, and find the elastic deformations and hierarchical 
atomic relaxations in the vicinity of a domain wall between two different homogeneous strain states. 
We clarify the microscopic origin of gradient terms, some of which are included phenomenologically 
in Ginzburg-Landau theory, by showing that they are anisotropic. 
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An understanding of atomic scale lattice distortions 
is essential for correctly describing the elastic energies 
of nano- structured materials. New generations of experi- 
mental tools to probe individual atoms and local environ- 
ments [lj , and the growing interest in complex functional 
materials, in which local lattice distortions are coupled 
to electronic, magnetic, and chemical degrees of freedom, 
further emphasize the need for a consistent theoretical 
framework to describe strain-based materials. For exam- 
ple, in perovskite manganites the change in oxygen ion 
displacement at each site is associated with the charge 
and orbital ordering states 2]. An atomic scale descrip- 
tion of the interface or domain wall between two different 
homogeneous states is thus a first step towards predicting 
functionality located at the domain walls. 

Strain variables (rather than displacement) with con- 
straints have been recently shown to have advantages 
for describing the long wavelength lattice distortions ob- 
served in, for example, martensitic materials and, more 
generally, solid-solid phase transformations The 
anisotropic long-range interaction in the order parame- 
ter strain fields drives the formation of a rich landscape 
of multiscale elastic textures. The aim of this work is 
to formulate a microscopic description of elasticity and 
demonstrate the relationship with and differences from 
long-wavelength continuum theory. We introduce ap- 
propriate inter-cell and intra-cell distortion modes and 
show how the form of the elastic energy recovers the 
correct phonon spectra. The discreteness of the lattice, 
choice of modes and constraints among them give rise to 
an anisotropic gradient expansion for the elastic energy. 
This leads to elastic domain wall solutions that are dif- 
ferent from those predicted using continuum theory; we 
obtain 0° and 90° 'staircase' domain walls for sufficiently 
small bulk modulus (or 'soft') materials, in addition to 
the 45° or 135° walls predicted from continuum theory 
for 'hard' materials. 

Models based on displacement variables with pair- 
potentials, such as Born- von Karman models |4|, have 
been widely used to incorporate 'microscopic elasticity'. 
However, essentially because distortion implies strain, 
the physical insight for atomic scale elasticity will reveal 



itself in the language of strain-related variables presented 
here. Moreover, our work is quite distinct from recent ef- 
forts to describe elasticity of nanometer-sized objects 0. 
The interest there is to describe long wavelength strains 
in a given dimension with other dimensions maintained 
at nanoscales, such as ultra-thin long nanowires. Our em- 
phasis is to describe atomic scale distortions, irrespective 
of whether the region of interest is in bulk or nano-sized 
objects. Our approach describes elastic deformation in 
terms of intra-ce\\ modes or "shuffles" of atoms, which 
are essential in describing short wavelength lattice dis- 
tortions, and distortion of unit cells, instead of adopting 
coarse graining approximations 0. 

We illustrate our ideas in detail for the simplest case, 
namely a square lattice in two-dimensional (2D) space 
with a monatomic basis. We find that the most conve- 
nient strain-related variables for atomic scale distortions 
are the normal distortion modes (more precisely, symme- 
try coordinates) of an elementary square object of four 
atoms, as shown in Fig. Because of the number of 
atoms in this object and the dimensionality, eight normal 
modes exist. The rigid rotation and the two rigid trans- 
lations (not shown in Fig. ^) cost no elastic energy, and 
therefore, are not distortion modes. The first three dis- 
tortion modes in Fig.[l]correspond to the usual dilatation 
(ei), shear (e2), and deviatoric (es) strains of the contin- 
uum elasticity theory for a square lattice 3]. The next 
two degenerate modes in Fig. ^ s + an d correspond 
to the "intracell" or "shuffle" modes of the square lat- 
tice [jj, which are absent in continuum elasticity theory. 
Our approach uses these five distortion variables defined 
for each plaquette of four atoms at z, z-h(lO), z-h(ll), and 
i + (01), where z represents the coordinate of the lattice 
points, to describe the elastic energy 8]. 

Since the five variables are derived from two displace- 
ment variables for each lattice site, they are related by 
three constraint equations. By representing e±, ei, es, 
s+, and S- in terms of displacement variables d x and d y 
in k (wavevector) space and eliminating d x and d y the 
constraint equations are obtained. One of them is the 
microscopic elastic compatibility equation, which relates 
strain modes: 
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FIG. 1: Normal distortion modes for a square object of four 
atoms in 2D. 



(1 — cos k x cos k y )e\(k) — sin k x sin k y e 2 (k) 
+ (cos k x — cos ky)es(k) = 0. (1) 
The other two relate the intracell and the strain modes: 
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±isin(^f^)e 3 (£) =0. (2) 
These constraints generate anisotropic interactions (re- 
flecting the lattice symmetry) between atomic scale 
strain fields, similar to the compatibility equations in 
Ref. 3], but now including the intracell modes. In the 
long wavelength limit, our description approaches the 
continuum model: For k — > 0, the above constraint equa- 
tions can be written in real space as 

V 2 ei (r) - 2V x V v e 2 {r) + (V* - V 2 x )e 3 (f) = 0, (3) 

s ± (f) = [(V„ ± V.) ei (r) + (V„ T V.) e 3 (f}} /4. (4) 
Equation (3) is the usual compatibility equation in con- 
tinuum theory. Equation (4) shows that the spatial varia- 
tions of strains always generate intracell modes, the mag- 
nitude of which vanish as the inverse of the length scale 
of the strain mode variations. It is well-known in con- 
tinuum Ginzburg-Landau theory that the energy associ- 
ated with the gradient of strains is responsible for do- 
main wall energies as, for example, in structural phase 
transitions Q. The above result shows that the intra- 
cell modes are at the origin of such energy terms. Since 
our strain-related variables become identical to conven- 
tional strain variables in the long wavelength limit, var- 
ious length scale lattice distortions may be described in 
a single theoretical framework. This makes it possible 
to study typical multiscale situations where both short- 
and long- wavelength distortions are important. It also 
provides a natural framework for incorporating interac- 
tions between atomic scale strain-related fields coupled 
to other degrees of freedom in functional materials. 

The following analysis of the simple harmonic elastic 
energy for the square lattice further exemplifies the util- 
ity of these variables. We consider the simplest energy 
expression by approximating the total elastic energy by 
the sum of the elastic energy of each square: 



n=l,2,3 



ra=+, 



where A n and B denote elastic moduli and 'intracell mod- 
ulus', respectively. The couplings between e±, e 2l es, s+, 
and S- at the same site are forbidden by symmetry at 
the harmonic level, but are allowed at the anharmonic 
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FIG. 2: An example of calculated phonon spectra (a) with 
and (b) without intracell modes for a 2D square lattice with 
a mono- atomic basis. The upper phonon branch is shown for 
both cases (M = h = 1). 



level, which may have important consequences for struc- 
tural "phase transitions" at the nanoscale. In Eq. (5) the 
inter-atomic elastic energies between atoms beyond each 
square, or further than the second nearest neighbors, are 
neglected. These interactions may be included by adding 
energy terms with distortion variables at different sites, 
e.g., ei(i)ei(i + (10)). Since some of the atomic pairs are 
shared by two square plaquettes of atoms, the parameters 
in Eq. (5) should be appropriately renormalized. A ro- 
bust way to determine the parameters is to compare the 
phonon spectrum of our model with experimental data. 

For the lattice energy of Eq. (5), the phonon spectrum 
is given by \fMhuj = yE^± \^E^, where E\ = (A\ + 
A2 + A$)(l — cos k x cos k y )/2 + B(l — cos k x )(l — cos k y ), 
E 2 = (A 1 + A 2 - A 3 ) 2 sin 2 k x sin 2 fc y /4 + (A 1 - A 2 + 
As) 2 (cosk x — cos/c 2y ) 2 /4, and M is the mass of an atom. 
A typical spectrum (upper branch) for Ai=5, A 2 =4, 
^3=3, and B=5 is shown in Fig. Efa). At k = (7r,7r), 
the distortion is a pure intracell mode, and the energy 
depends only on the intracell mode modulus B. There- 
fore, as shown in Fig. Hfb), cj(7r,7r) vanishes without the 
intracell mode (5=0), which is unphysical. As k — > 0, 
the slope of the phonon spectrum is determined only by 
strain mode moduli, A±, A 2 , and A3, whereas the cur- 
vature depends on B as well because of Eq. (4). For 
k — > 0, since the intracell modes vanish as the in- 
verse of the wavelength, the lattice energy approaches 
■E sq .iat(fc -> 0) « EfeVi,2,3^ e nW/ 2 5 in agreement 
with continuum theory. Using Eq. (4), the energy for 
the intracell modes in Eq. (5) can be written as i^intra ~ 
/ drBKVe^ 2 + (Ve 3 ) 2 + 2(W y e 1 V y e 3 - V x e x V*e 3 )]/16. 

We compare our approach to a displacement-based 
Born- von Karman model |^ for the square lattice. The 
first-nearest and the second-nearest neighbor atoms are 
connected by central-force and non-central-force springs 
with spring constants iff, iff, iff and iff, re- 
spectively. Elastic energies of the lattice, one from 
our model and the other from the Born-von Karman 
model, can be represented in the following form in terms 
of displacement variables, d x (i) and d y (i): E sq .\ at = 
Ek, a ,b da (-k)Dab(k)d b (k). We find that the £> a6 (fc)'s 



for both models become identical if: A\ = iff 
2{K% - K* ), A 2 = 2(K( 



N 



K?), A 3 = K 



c 



3 



K? + AKg, and B = K° + iff. 

We apply our formalism to obtain the domain wall so- 
lution for the atomic displacements between two homo- 
geneous strain states (a "twin boundary" ) due to a phase 
transition to a rectangular lattice. We then compare the 
solution to that obtained from continuum theory where 
discreteness effects are neglected Q • With elastic energy 



E^L where 



= Y,i\Mei$) 2 + \A 2 e 2 $) 2 + \B[s + (i) 2 +S-$)% 

E%1 = J2,-lA' 3 e 3 (i) 2 + ±F 3 e 3 (i) 4 , (6) 
the degenerate ground state of E rec is a uniform state 
with e3=± y/A' 3 /Fs, and ei=e 2 =s+=S-=0. To study 
the domain wall between these two degenerate rectan- 
gular ground states, we consider es(i) as the order pa- 
rameter and minimize Erll with respect to the other 
variables using the constraint equations [Eqs. (1) and 
(2)] and the method of Lagrange multipliers. We obtain 
£ r ( ec,min = E fc - 5 e 3 (~k) U(k)e 3 (k) , where U(k) = (V, + 
V 2 + Vs)/Vi, and V 1 = [A 2 (A 1 /3i + A 2 p 2 ) + B 2 0l)Ai(3 3 , 
V 2 = (2A 1 A 2 p 2 + A 2 0l + A 2 (3 2 )B^, V 3 = 2{A l (3 2 + 
A 2 0i)B 2 0\ + B 3 p 1 pl Vi = (A 1 (3 2 + A 2 (3 2 + Bfcfa) 2 , 
with fa = 1 — cos k x cos k y , fa = — sin k x sin k y , fa 
cos k x — cos k y , and fa = (1 — cos fc £C )(l — cos fc^). 

With k x =k cos 9 and k y =ksm6, the expansion of 
f7(fc,0) about fc=0 yields U(k,0) = U o {0) + /7 2 (<9)/c 2 + 
O(V), where f7 o (0) = A X A 2 cos 2 20/ (Ai sin 2 2(9 + A 2 ), 
and /7 2 (<9) = sin 2 2(9 [6 AiA 2 5 sin 2 2(9 + 4AiA 2 (Ai + 
A 2 ) cos 2 20 + 3B(A| + A 2 sin 2 20)]/[24(A 2 + A sin 2 2<9) 2 ]. 
The term C/ is purely orientation-dependent without a 
length scale, and is minimized at 6 =45° and 135°, as ob- 
tained in Ref. 3] . The difference between continuum and 
our discrete theory lies in the k 2 term: continuum the- 
ory commonly assumes isotropic gradients in the order 
parameter, i.e., (Ves) 2 7], whereas U 2 (Q) is anisotropic. 
The two origins of the anisotropy are: (a) the compati- 
bility relation, Eq. (1), which has higher powers in k than 
Eq. (3) due to discreteness, and (b) the presence of shuf- 
fle mode energy. The latter can be written as gradients 
of strains, but with corrections to the phenomenological 
isotropic term, (Ves) 2 , used in Ginzburg-Landau theory. 
As U 2 {6) is minimized for = 0° and 90°, it competes 
with U o (0) which prefers = 45° and 135°. Thus, the 
domain wall direction depends on the length scale with 
a critical length scale A c ~ y/B/Ai. If A c < 1, i.e., less 
than the interatomic spacing, the domain wall has direc- 
tion 45° or 135° down to atomic scales. If A c > 1, then 
for length scales smaller (larger) than A c , the domain wall 
direction is 0° or 90° (45° or 135°) and the domain wall 
has multiscale attributes. 

We examine first the case A c < 1 that would ap- 
ply to materials with relatively large bulk modulus A\ 
('hard' materials) for fixed B. Here k x = ±k y and 
U(k) = B(l — cosfc ;E )/(l + cosk x ). The domain wall 
width is a result of the competition between U(k) that 
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FIG. 3: (Color) Atomic scale 135° domain wall profile for 
critical length scale, A c < 1 along the direction perpendic- 
ular to the domain wall: (a) strain e 3 and shuffle s_, (b) 
differences in e 3 (Se 3 =e 3 , atomic - e 3 , continuum), s- (Ss-) and 
displacement parallel to the domain wall direction (Sd\\) be- 
tween the results from continuum theory for k ~ and our 
model that includes discreteness. The fields ei, e2, s+, and 
displacement perpendicular to the domain wall are zero. Pa- 
rameter values are A\ — 5, A2 = 4, A' 3 = 4, B = 5, and 
F 3 50. 
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favors k — > 0, or thick domains, and E^c that favors 
sharp walls. We illustrate the domain wall solution with 
135° domain wall direction. The only non-zero distor- 
tion modes are es and s_ (s+ for a 45° domain wall). 
The strain es reverses sign at the domain wall, the in- 
tracell mode S- is confined within the domain wall, and 
the atomic displacements are parallel to the domain wall 
direction. The numerical solution for es and S- along 
a line perpendicular to the wall is shown in Fig.|3fa), for 
which A c ~ 1. (Narrow domain walls with widths of a 
few unit cells, as considered here, have been identified ex- 
perimentally 10].) The corresponding displacement field 
near the center of the domain wall is shown in Fig.0fa), 
in which the red and blue colors show regions with es pos- 
itive and negative, respectively. Both figures show that 
the center of the domain wall is located at bonds rather 
than sites to avoid the higher energy state of es=0 and 
large s_. As for a Peierls-Nabarro barrier the higher 
energy (by 4.4 x 10 -4 per unit length for our parameter 
values) for the site-centered domain wall acts as a pinning 
potential for the domain wall due to the inherent discrete- 
ness. In Fig. 3(b) we compare our results with contin- 
uum theory, which predicts es = e™ 3 "* tanh(z s /£) [3 and 
8_ = des/2di s from Eq. (4), where i s = i x + i y . The dif- 
ferences in the interface region, shown in Fig.|3fb), are of 
the order of 10% of ef a * = y/A!jF$. The domain wall 
width |lOj is roughly given by 2£ = V2y/B/A' 3 and the 
ratio between the maxima of s_ and e3, s^^/e^, is 



about l/(2f). 

Anisotropic effects in U(k,0) become more apparent 
away from equilibrium, e.g., at finite temperatures or 
in other conditions where met ast ability is present. Fig- 
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FIG. 4: (Color) Atomic displacements in the vicinity of do- 
main wall: (a) A c < 1 and (b) A c > 1. Color represents 
the sign of ez (red: positive, blue: negative), and is lighter 
compared to Figs. |5] and [SI to make the arrows visible, (c) 
The displacement component perpendicular to the large scale 
domain wall direction for (b), magnified by a factor of 3. Pa- 
rameter values are A2 = 4, A' 3 = 4, B = 5, F3 = 50, and A\ 
=5 for (a) and A±=l for (b) and (c). 
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FIG. 5: (Color) Nonequilibrium domain wall state: (a) 33, 
(b) s+, (c) s- for A c < 1, (d) ez for A c > 1. Parameter values 
are identical to Fig.0] Dark red corresponds to 0.28 and dark 
blue to -0.28. Green implies a value close to zero. 



ure EJa) shows the results for e% of a 2D simulation 
away from equilibrium , in which atomic scale domain 
walls are oriented along 45° and 135° directions. The 
corresponding intracell modes are shown in Figs. Etb) 
and[H{c). Note that both s + and s_ shuffles are present 
at or near interfaces and only one of these modes sur- 
vives at equilibrium. The horizontal or vertical "jogs" 
in the s+ and S- walls are secondary defects due to the 
competition between Uq and L^, which provide principal 
relaxation forces. 

The domain wall solution for A c > 1, typical for small 
bulk modulus A\ or 'soft' materials, is shown in Fig. 
for which A c ~ y/b. The e% field in Fig. EJa) shows that 
on length scales of the size of the system (larger than 
A c ), the diagonal orientation is still preferred. However, 
this diagonal domain wall consists of a 'staircase' of 0° 
and 90° domain walls of length scale A c . The existence of 
0° and 90° walls in e% forces elastic compatibility to in- 
duce alternately large positive and negative values in the 
dilatation strain e\ in the horizontal and vertical parts 
of the 'staircase', as shown in Fig. EJb). This has im- 
plications for the functionality of the domain walls. For 
example, the e\ field can couple to charge and modulate 
the local charge density along the wall. Similar features 
are also reflected in s + and s_, as shown in Figs. Efc) 
andEJd). The displacement pattern within the square in 
Fig. |H|is shown in Fig. Hfb). Unlike the case A c < 1 in 
Fig. ^a), the displacement has a component perpendic- 



FIG. 6: (Color) Atomic scale domain wall solution for materi- 
als with A c > 1. Parameters are the same as in Figs.lUb) and 
H{c), which show the region inside the square in this figure. 
Strain e2 is zero. Color scheme is the same as in Fig. 

ular to the 135° large scale domain wall direction, which 
is shown in Fig. 2fc) magnified by a factor of 3. This 
component is greatest for the atoms at the boundary of 
the two domains. A nonequilibrium state is shown in 
Fig-Efd). The small domains have 0° or 90° boundaries, 
but over larger length scales these domains are correlated 
along 45° and 135° directions. 

In summary, we have reported an approach to "atomic- 
scale elasticity" , which uses symmetry modes of elemen- 
tary objects of atoms as distortion variables. A gradient 
expansion for the energy with anisotropic coefficients has 
been obtained, with corrections to the usual phenomeno- 
logical isotropic gradient terms used in Ginzburg-Landau 
theory. As an illustration, we have obtained domain wall 
(twin boundary) solutions in terms of strain and intra- 
cell modes and have shown how they differ from the con- 
tinuum elastic soliton solution 7]. Our work provides 
the basis for interpreting atomic scale features in HREM 
studies of domain walls . 
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